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ABSTRACT 

I examine the angular momentum evolution during the 1837-1856 Great Erup- 
tion of the massive star rj Carinae. I find that the new estimate of the mass blown 
during that eruption imphes that the envelope of t] Car substantially spun-down 
during the 20 years eruption. Single-star models, most of which require the en- 
velope to rotate close to the break-up velocity, cannot account for the bipolar 
nebula— the Homunculus— formed from matter expelled in that eruption. The 
kinetic energy and momentum of the Homunculus further constrains single-star 
models. I discuss how rj Car can fit into a unified model for the formation of 
bipolar lobes where two oppositely ejected jets infiate two lobes (or bubbles). 
These jets are blown by an accretion disk, which requires stellar companions in 
the case of bipolar nebulae around stellar objects. 

Key words: binaries: close— circumstellar matter— stars: individual: tj Carinae— stars: 
mass loss— stars: winds 



Introduction 



The bipolar structure of the Eta Carinae [t] Car) nebula— the Homunculus— is not 
unique. Its basic structure, that of bubble pair, where each bubble is bounded by a thin 
dense shell, filled with a low density interior, and having a narrow waist between them, is 
seen in many diverse objects. Examples include the Perseus cluster of galaxies (Fabian et 
al. 2002), the symbiotic nebula He 2-104 (Corradi & Schwarz 1995; Corradi et al. 2001; here 
the dense shell is not closed), and the planetary nebula (FN) NGC 3587 (FN G148.4+57.0; 
e.g., Guerrero et al. 2003). The similarity between bubble pairs in clusters of galaxies and 
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in PNs is discussed in Soker (2003a,b; 2004a), while in Soker (2004b) I discuss the similarity 
of these systems with symbiotic nebulae and the bipolar nebula of rj Car. In clusters of 
galaxies such bubbles are known to be formed by oppositely ejected jets, which are detected 
by radio emission (e.g., Hydra A; McNamara ct a. 2000). In a series of papers I argue that 
this suggests that bubble pairs in PNs (Soker 2003a,b; 2004a) and symbiotic nebulae (Soker 
2004b) are also formed by such jets (note that not all PNs are shaped by jets). In partic- 
ular, there are more and more observations showing and hinting at the presence of jets in 
symbiotic systems (e.g., Kellogg, Pedelty, & Lyon 2001; Brocksopp et al. 2004), in PNs, and 
PN progenitors (e.g., Imai et al. 2002; Hirano et al. 2004; Sahai et al. 2003). 

In previous papers I attribute the bipolar structure of 77 Car (e.g., Ishibaslii et al. 2003) 
and the presently blown fast polar wind found by Smith et al. (2003a) to interaction with 
the binary companion (Soker 2001, 2003c). In particular, in Soker (2001) I proposed that 
the bipolar nebula was shaped by jets blown by the companion, via an accretion disk, during 
the 20 years Great Eruption a century and a half ago. In that paper I listed some arguments 
in favor of such a model. However, the new finding of Smith et al. (2003b) of a more massive 
Homunculus, and several new papers using a single star model for the shaping of the wind and 
circumstellar matter of f] Car (e.g., Dwarkadas & Owocki 2002; Smith et al. 2003a; Gonzalez 
et al. 2004; van Boekel et al. 2003) motivate me to reconsider the single star model. Other 
papers in recent years study the role of rotation and/or binary companion in 77 Car. Stothers 
(1999) find that the rotation does not affect much the instability of luminous blue variables 
(LBV), which is in support of the view I take in the present paper. Davidson (1999; see 
also Davidson 2000), on the other hand, argues that a binary model for 77 Car is not "a 
Panacea". I accept Davidson's view that some crucial aspects of the r) Car eruptions arc due 
solely to the behavior of the massive star, and not to the companion. My main interest is in 
showing that the bipolar structure of the nebula is due to a binary companion. In section 
2 I study the momentum and energy budget of f] Car. In section 3 I consider the required 
stellar angular momentum in single-star models. I will not follow the evolution of rotating 
massive stars with any specific model (for these see, e.g., Heger & Langer 2000; Maeder, 
& Meynet 2003; Meynet & Maeder 2003). I will use analytical approach. In section 4 I 
summarize my finding that single-star models encounter severe problems. 



2. Wind's energy and momentum 

Some of the arguments given here are not new. However, I put them in a broader context 
in light of the new estimated mass that was blown during the 20 years Great Eruption 
Mge — I2M0 (Smith et al. 2003b). For a wind driven by radiation pressure the effective 
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number of times a photon is scattered by the wind material, i.e., in the positive radial 
direction, is given by 

^M^^ / u yV ^ ^ ( ^ m 

~ L/c \0.5Mq yr-i j UoO km s-V ' ^ ' 

where My, is the mass loss rate into the wind, v^, is the terminal wind speed, is the 
stellar luminosity, and c the light speed. The different variables are scaled with the mass 
loss rate (Smith et al. 2003b) and luminosity (Davidson & Humphreys 1997) during the 
Great Eruption. This is a very large number, which is not encountered in other stars. For 
example, in most PNs and highly evolved red giant stars < 1 (Knapp 1986). The PN 
with the highest value of rjg ~ 600 in the list given by Kanpp (1986) is NGC 2346. The 
central star of this PN contains a close binary companion (Bond 2000), which went through 
a common envelope evolution and released large amount or orbital energy. Therefore, the 
binary companion is behind the large value of rfs in the PN NGC 2346. Smith et al. (2003b) 
notice that such a high value of r]s occurs in explosions. However, the mass loss rate during 
the Great Eruptions lasted for ~ 20 years, a time longer by more than an order of magnitude 
than the dynamical time of r] Car, even if it swelled to ~ 10 AU. Shaviv (2001) compares the 
super-Eddington wind of the Great Eruption to that of novae and obtains satisfactory fit to 
the Great Eruption outflow, assuming spherical wind. Although he assumes a steady state 
mass loss process, the photosphcric radius of the nova ejected mass expands by several orders 
of magnitude during the explosion, and the mass involved occupies a thin layer on the white 
dwarf surface prior to the explosion. The ~ lOM© ejected mass in the Great Eruption, on 
the other hand, originated in a thick envelope layer of r] Car. In addition, the formation of 
the bipolar structure in an explosive event requires different models than the one I criticize 
in the present paper. The results of Shaviv (2001), though, show that a single star can lose 
mass at a high rate. Such a high mass loss rate is needed in the presently proposed binary 
model, but here the wind speed, hence the momentum and kinetic energy supplied by the 
primary, can be much lower. 

The new estimated energy of the Homunculus is -Enom ~ 10*^^-^ — 10^° erg (Smith et al. 
2003b). The ratio of this energy to the energy radiated during the 20 years Great Eruption 
is ^ 

_ Enom _ 1 ( \ I Enoia 

^ X 20 yr ~ VIO^-^Lq J \5 x lO^s erg 

This imphes that half the energy hberated during the Great Eruption was radiated, and 
half was converted to the kinetic energy of the wind. The present mass loss rate of < 
2 X 1O~^M0 yr~^ (van Boekel et al. 2003) over the ~ 150 years span since the Great Eruption 
contributes negligible amount to this energy. 
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The required transfer efficiency of momentum and energy from radiation to the wind in 
single star models of rj Car is much higher than in any other object blowing similar winds. 
These values are more typical for binary systems, as I now discuss. In principle there are 
two processes by which a binary companion can account for such a high kinetic energy of the 
wind. First, the binary system can release orbital energy AEoa — 0.5GMiM2/af where Mi 
and M2 are the two masses inward to the final orbital separation of a/, and assuming that 
the final orbital separation is much smaller than the initial orbital separation. If a fraction 
X of the released orbital energy is deposited into the expelled mass, then the final orbital 
separation for a companion to explain the Homunculus energy is 

- (5 x^Perg) (^) (12^) (3^) 

The presently observed orbital period is 5.5 yr, which for primary and companion masses 
of Ml = I2OM0 and M2 = 30Mq, respectively, implies an average orbital separation of 
a ~ 16.5 AU (for observational support for the presence of a binary and its properties see, 
e.g., Damineli 1996; Ishibashi et al. 1999; Damineli et al. 2000; Corcoran et al. 2001a,b; 
Pittard & Corcoran 2002; Duncan & White 2003 Fernandez Lajus et al. 2003). Therefore, 
the orbital energy of the presently observed binary system cannot account for the kinetic 
energy of the Homunculus. However, 77 Car could have swallowed a closer, lower mass third 
star during its Great Eruption. For example, a third star with M3 = 5Mq could have spiral 
down to the core of rj Car, and released more than the required energy. The Lesser Eruption 
of 1890 makes this scenario unlikely, although it can't be ruled out based on pure physical 
arguments. 

In the second process, the companion is outside the primary envelope, it accretes mass 
via an accretion disk and blows two jets (or a collimated fast wind: CFW). In this scenario, 
which was proposed for the formation of the Homunculus in Soker (2001; where more sup- 
porting arguments are given), the primary in rj Car expelled its mass at low speeds. In 
principle, the companion may accretes half of the mass expelled by the primary and blow a 
fraction of ~ 0.2 at a speed equals to its escape velocity Vj ~ 2000 km s~^. For an accreted 
mass of ~ IOMq the energy carried in the jets is then ~ 8 x 10^^ erg, enough to account 
for the energy of the Homunculus. In addition, in such a scenario a substantial fraction of 
the energy radiated by the system during the great eruption came from the accretion energy 
onto the companion. Recently, Smith & Morse (2004) reported the discovery of extremely 
fast material, v > 3200 km s~^, in 77 Car. They attribute this material to mass loss from the 
primary star. In the binary model, on the other hand, this extremely fast material was first 
accreted onto an accretion disk around the companion, and then was ejected at this high 
speed. 
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3. The angulcir momentum evolution of rj Car 

3.1. Moment of Inertia 

Following Soker & Harpaz (1999) ^ I examine the ratio of the density at the photosphere, 
Pp, to the average envelope density Pa- As shown below, this ratio indicates the moment of 
inertia, a relevant quantity for the slowing down process. The photospheric density is given 
by (Kippenhahn & Weigert 1990) 

3 kB R-'nr ^ ^ 

where iiuih is the mean mass per particle, /c^ is the Boltzmann constant, and k is the opacity. 
Substituting typical values for 77 Car in the last equation gives 

where a black body luminosity was assuemd to eliminate the stellar radius. Using the (solar 
composition) opacity as given by Alexander & Ferguson 1994; see also Rogers & Iglesias 
1992), I find the following approximate, but adequate, fitting in the relevant photospheric 
densities 

_ r 10(r/10^ K)^^ cm^ g-i for 5, OOOX < T < 8, 500X 
'^~\ 1 cm2 g-i for 8500X < T. ^ ' 

I take the envelope mass M^^y to contain about half of the 77 Car stellar mass. For 
LBV as 7] Car, the envelope mass is expected to be lower, hence increasing the slowing down 
process and strengthening the effect studied here. The average density in the envelope is 
given by 

As was shown for AGE stars (Soker & Harpaz 1999), and can be checked for more massive 
stars, when pp is not much smaller than pa, the envelope density profile is shallow. In low 
mass envelopes, the density profile is almost fiat in the outer region (Soker & Harpaz 1999). 
For massive stars this can be seen in the model of an initial I2OM0 star which was reduced 



^Note that the density scale in Figs. 1-5 of SH99 is too low by a factor of 10; the correct scale is displayed 
in their Fig. 6. 



-6- 



to 66.6M0 as calculated by Stothers & Chin (1993). In that model (their fig. 2), the effective 
temperature is Tp = 8000 K, and luminosity 2 x IO^Lq. The radius is ~ 700Rq. The envelope 
mass is very low, because of the mass lost by the star. On average, the density profile from 
~ 0.05-R^, to the surface is p oc r~^. The parameters in the evolved 120Mq model of Maeder 
(1981) without mass loss are L ~ W^Lq, Tp ~ 4000 K, i? ~ 20001?©, and Mc„v = SOM©. 
From the equations above, and using the correct opacity for low densities of k ~ 3 x 10~^ 
I find pp — 10~^° g cm~^ and pa = 5 x 10~^ g cm~^. The ratio of 50 is moderate, and the 
density profile in the outer envelope is Pe ^ 

The structure of the envelope considered above determines the moment of inertia of the 
star (that of the core is negligible for giants) 

/ = aMenvi?^. (8) 

For an envelope density profile of Pe oc one finds a — 2/9, while for pe oc the value 

is a = [3 ln(i?/rin)]~^, where rin is the inner radius of the envelope. In the model described 
above of Maeder (1981), -R/nn ~ 10, and I find a ~ 0.15. For steeper density profiles 
the value of a is lower. In the sun, for example, equations (5)- (7) give a very large ratio 
of log(pa/pp) ~ 5. The density profile in most of the solar interior (beside the very outer 
envelope) can be fitted by 

p(r) = poCxp{-Kr), (9) 

with KRq= 10.54 (Bahcall, Pinsonneault, & Basu 2001). For such a profile, and with 
Tin <^ -R, the moment of inertia coefficient is a ~ 8/{KR)'^. For the sun this gives a = 0.07. 

This subsection shows that the ratio p^/ pp can be used as a crude indicator for the 
moment of inertia. The exact value of a can't be predicted from this ratio, but for the 
present goal it is enough to use the crude relation 

r 0.2, for log(p„/pp) < 2 
ar^l 0.1-0.2 for 2 < log(p„/pp) < 4 (10) 
[ <0.1 for4<log(p„/p,). 



3.2. Angular Momentum Loss 

Consider an envelope of a giant star rotating as a solid body. The angular momentum 
loss rate from the envelope to the wind is 

jwind = PUR'M, (11) 

where cu, is the stellar angular velocity, J is the stellar angular momentum, and /3 depends on 
the mass loss geometry: for a constant mass loss rate per unit area on the surface (3 — 2/3, 
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while for an equatorial mass loss P — 1. Smith et al. (2003b) argue for an enhanced polar 
mass loss rate during the Great Eruption, for which /5 < 2/3. However, still some extra 
mass resides in the equatorial plane even according to Smith et al. (2003b). Therefore, I will 
take f3 ~ 0.5 — 0.7 for the Great Eruption mass loss geometry. The change of the envelope's 
angular momentum with mass loss is given by (e.g., Soker & Harpaz 2000) 

dlnJenv dlnuj dlnl /3 

+ T^ — ^7 — ~ ~TT7 — T = \^^) 



dlnMenv cilnMenv din M^nv Q;(Menv) 

If the structure of the atmosphere does not change much while mass loss occurs then 
dlnl/dlnMenv — 1 and 5 is constant. The solution of the last equation becomes 



CO /M„ 



'env 



'envO 



(13) 



3.3. Anguleir Velocity Evolution During the Great Eruption 

Most single-star models I am aware of for the formation of the Homunculus require r] 
Car to rotate at Q = (cj/cjkcp) ^ 0.7, where uj is the angular velocity of the stellar envelope 
and ci^Kep is the Keplerian velocity on the equator (i.e., the break-up angular velocity). 
Dwarkadas & Owocki (2002) and Smith et al. (2003a), for example, take Q — 0.9, while 
Maeder & Desjacques (2001) take Q = 0.8 — 0.9. The model presented by Langer, Garcia- 
Segura, & Mac Low (1999) is different in that they consider the ratio of luminosity to the 
Eddington limit. Similar considerations are qualitatively presented by Zethson et al. (1999) 
to account for the slow ejecta in the equatorial plane of rj Car. They emphasize that even 
slowly rotating luminous stars can have highly non-spherical mass loss geometry. They 
also point out that their model cannot explain why the slow ejecta are in dense compact 
condensations. In the binary model, these condensations are assumed to be formed by the 
influence of the companion. The model of Langer et al. (1999) which gives two lobes similar 
to those in rj Car have Q < 0.2. However, most of the mass in their model is being lost 
at low velocities, < 400 km s~^, and it resides in the equatorial plane rather than in the 
Homunculus. Therefore, the total kinetic energy in their model is much below that of the 
Homunculus. Aerts, Lamers, & Molenberghs (2004) study the influence of rotation on the 
mass loss from r] Car, and also find the wind velocity to be too low. Therefore, in what 
follows I will refer only to the models require Q >0.7. 

To facilitated an analytical treatment that will demonstrate the problems of a rotating 
single star, I consider two extreme cases: an eruption accompanied by expansion of the 
envelope to i? ~ 10 AU (e.g., Davidson & Humphreys 1997), and an eruption leaving the r) 
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Car primary star a hot star. The angular momentum of r) Car at the beginning of the Great 
Eruption is determined by both single star evolution and binary interaction. To posses fast 
rotation, though, the star must have interacted with a stellar companion, because single stars 
slow down with mass loss (Maeder & Meynet 2000), e.g., large mass loss events before the 
Great Eruption, such as the one proposed by Bohigas et al. (2000), would have substantially 
slowed down r] Car. As a star like rj Car expands to ~ 10 AU, its moment of inertia 
increases both because of the increase in size and the increase in moment of the inertia 
coefficient a (eq. 8). Presently, rj Car has a radius of ~ lOOi?© and an effective temperature 
of ~ 20, 000 - 30, 000 K (Hillier et al. 2001; Smith et al 2003a). For such a model a < 0.1 
by equations (5), (7) and (10). For a Great Eruption luminosity of lO^-^L© and a radius 
of ~ 10 AU, the effective temperature becomes ~ 8,300 K, and I find a ~ 0.1. Such an 
expansion implies, because of conservation of angular momentum, that the envelope must 
rotate very slowly, even if before expnation the star was almost at break-up angular velocity. 
I take now the following parameters: total mass lost during the Great Eruption equal to 
AMge = I2M0; an envelope mass of about half the stellar mass, or Mgnvo ~ 70Mq at the 
beginning of the great eruption (as stated above, this is an upper limit on the envelope 
mass; the envelope mass is likely to be lower, increasing the effect studied here); mass loss 
geometry with /3 = 0.5 (see discussion following eq. 11); and a = 0.1 — 0.15. For these 
parameters, S — 1 = {f3 / a) — 1 = 4. By equation (13), the envelope angular velocity at the 
end of the Great Eruption is f2/ = i^o (58/70)^ = 0.5^0, where Qq is the (non-dimensional) 
angular velocity before mass loss starts but after expansion. At the middle of the Great 
Eruption, Q^n = ^^0(64/70)"^ = 0.7Qo. For a = 0.15, as in the extended 60Mq model of 
Maeder (1981), these values are flj = 0.65r2o and flm = O.Sflo. Since in this scenario the 
angular velocity at the beginning of the Great Eruption should be very low as well ^ 1, 
I conclude that in the case when 77 Car expanded during the Great Eruption the angular 
velocity is very low, and a single star model cannot account for the bipolar structure. 

Consider now the other extreme, where rj Car remained a hot star during the Great 
Eruption. Equations (5) (7) and (10) implies a < 0.1. Again, I take (3 = 0.5, i.e., a moderate 
polar enhanced mass loss rate. The value of P cannot be too low. This is because a low value 
of (3 implies highly enhanced mass loss rate along polar directions. This, however, will make 
the problem of the radiation momentum transfer to the wind (eq. 1), even more severe, as 
only radiation escaping along the polar directions accelerate the mass there. Therefore, a 
value of (5 = 4, as in the extended envelope case, may be an underestimate. In any case, 
for this value of 6, and for a maximum, possible value of f2o = 1, the angular velocity and 
middle and end of the Great Eruptions are flm = 0.7 and = 0.5, respectively 

The main conclusion of this section is that as a single star, rj Car could not have 
maintained a fast rotation during the 20 years Great Eruption. 
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4. Summciry 

Most single-star models for the formation of the bipolar nebula— the Homunculus— of rj 
Car are based on fast rotation of the progenitor during the 20 years Great Eruption, 1843- 
1856. These models do not consider the origin of the fast rotation, or the way the envelope 
maintains its fast rotation during the Great Eruption. A single-star spin-up mechanism as 
proposed by Heger & Langer (1998) requires the outer convective part of the envelope to 
loss mass to inner regions. This will decrease substantially the ratio Menv/Menvo in equation 
(13), hence substantially reducing the angular velocity as a result of the mass loss. The 
single-star models for the bipolar Great Eruption, therefore, do not directly address the 
question of what is the main factor behind the bipolar structure. There is, however, more 
to the single star model. Smith et al. (2003a), for example, attribute the 5.5 year period to 
angular momentum redistribution in the envelope. The star eruptively expels large amount 
of mass and angular momentum, such that its surface slows down. The star then relaxes, 
and angular momentum diffuses outward. After a period of 5.5 years, the surface rotates 
fast enough for a new eruption to occur. Although I see some problems with this scenario, it 
does, however, show that the angular momentum evolution in a single star scenario is more 
complicated than what I described in previous sections. 

The goal of the present paper was to show that the origin of the fast envelope rotation 
required by these models cannot be ignored. Even when fast rotation occurs, it is not clear 
how much it can affect instabilities and the mass loss process (Stothers 1999) I used general 
considerations, rather than a specific model. (Most papers dealing with a fast rotating 
progenitor of rj Car do not actually provide such a stellar model!) The first problem is the 
initial fast rotation required. This in principle can be provided by a relatively low mass 
companion which entered the envelope of t] Car in the past. The more difficult problem to 
overcome, as I showed in section 3 above, is to maintain a fast envelope rotation during the 
Great Eruption itself. On top of the angular momentum problem, the models are further 
constrained by the huge momentum and kinetic energy of the mass blown during the great 
eruption (section 2 above). These, for example, rule out a model where most of the mass is 
blown from the polar caps of the progenitor. 

Binary companions, on the other hand, can relatively easily account for the shaping 
of the Homunculus (Soker 2001), and the energy and momentum of the Great Eruption 
(section 2 bove). The basic process is that of a companion accreting from the primary 
wind, forming an accretion disk, and blowing two opposite jets. The energy source is the 
gravitational energy of the mass accreted on the companion. Comparing LBV's nebulae 
with other objects, e.g., PNs, strongly suggests that binary interaction is indeed behind the 
shaping of these non-spherical nebulae (e.g., O'Hara et al. 2003). The binary model has 
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another advantage: it incorporates rj Car into a unified model explaining all objects having 
two low density lobes with an equatorial waist between them; sometimes these are termed 
bubble pair. These objects include clusters of galaxies, where the bubbles are X-ray deficient 
bubbles, symbiotic nebulae, and planetary nebulae (see Soker 2003 a,b; 2004a). The shaping 
in all these objects is via supersonic jets. The jets are blown by accretion disks at a speed 
about equal to the escape velocity from the accreting object (Livio 2000). 

Following the analytical exploratory papers of the shaping problem from a binary-model 
perspective (Sokcr 2001; 2003c, and the present paper), the next step in understanding the 
shaping of the Homunculus is to conduct a 3D numerical simulation of the mass transfer from 
r] Car to its companion, including two oppositely ejected jets, with possibly wide opening 
angle (Soker 2004a), and including the orbital motion. 
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